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We have investigated the electromagnetic decays of the antitriplet and sextet charmed baryon 
systems with J p = | + , | + in the framework of the heavy baryon chiral perturbation theory. We 
first construct the chiral Lagrangians at 0(p 2 ) and 0(p 3 ). Then we calculate the electromagnetic 
(EM) decay amplitudes of the charmed baryon systems up to 0(p 3 ). With the help of the quark 
model, we estimate the low energy constants. The numerical results of the EM decay widths show 
good convergence of the chiral expansion. We notice that the two neutral EM decay processes 
S(,° —> 7 + E° and H* ,Q —> 7 + are strongly suppressed by the SU(3) U-spin flavor symmetry. 

With the same formalism, we also estimate the EM decay widths of the bottomed baryons. The EM 
decay widths of the heavy baryons may be measured at facilities such as LHCb and JPARC. The 
explicit chiral structures of the heavy baryon decay amplitudes derived in this work may be useful 
to the possible chiral extrapolations of the future lattice simulations of these EM decay amplitudes. 

PACS numbers: 12.39.Fe, 14.20.Lq, 13.30.-a, 13.40.Hq 


I. INTRODUCTION 

The low energy behavior of quantum chromodynamics (QCD) is highly nonperturbative and complicated. The 
study of the hadron properties helps us explore the infrared structure of QCD. In the past decades many heavy 
baryons were discovered experimentally |lH3| . 

The heavy flavored baryon systems containing a single heavy quark are particularly interesting. In the heavy quark 
limit, there exists the additional heavy quark spin and flavor symmetry when the heavy quark mass goes to infinity. 
The strong dynamics within the heavy baryon is governed by the two light quarks while the heavy quark acts as a 
color source and the spectator. The physical observables can be expanded in terms of 1 /toq where toq is the heavy 
quark mass. 

In this work, we focus on the heavy baryons containing a single charm or bottom quark and investigate their 
electromagnetic (EM) decays. The heavy baryon ground states contain the spin-i sextet, antitriplet and the spin-| 
sextet. For brevity, we refer to them as the 6, 3 and 6* dimensional baryons in the following. We will explore the 
property of the electromagnetic transitions between the heavy baryon multiplets. 

In literature, various theoretical approaches were employed to study the magnetic moments and EM decays of the 
light baryons [ ll 1 2(1. We have performed a systematic study of the chiral corrections to the heavy baryon mass and 
axial charge in Ref. [2j|. In this work, we will investigate the electromagnetic decay width of the antitriplet and 
sextet heavy baryon systems in the framework of heavy baryon chiral perturbation theory (HBChPT). We truncate 
the chiral expansion at 0(p 3 ). The tree-level electromagnetic decay amplitude appears at 0{p 2 ) and 0(p 3 ). We 
include the chiral loop corrections up to 0(p 3 ). The recoil corrections in terms of jj appear at 0(p 4 ) or even high 
order. As will be seen below, our results show good convergence of the chiral expansion. 

There are twelve possible EM decay channels. Because of the constraint of flavor structure, only ten channels 
arise at the leading order. At 0(p 3 ), there appear six additional EM interaction low energy constants involving the 
6-dimensional and 6*-dimensional baryons at the tree level. In contrast, the loop contributions appear at 0(p 3 ) and 
contribute to all twelve channels. 

The EM width is proportional to the cubic of the photon momentum, which is roughly equivalent to the mass 
difference of the initial and final states. The average value of the mass difference is 127 MeV between the 3 and 
6-dimensional charmed baryons, 67 MeV between the 6 and 6*-dimensional baryons, 194 MeV between 3 and 6*- 
dimensional baryons. In the calculation of the EM decay widths, we use the experimental values of the charmed 
baryon masses as input to derive the photon momentum. The EM decay widths of different channels differ greatly, 
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mainly due to the large difference of the photon momentum. In the appendix, we also extend the same formalism to 
estimate the EM widths of the bottomed baryons. 

This paper is organized as follows. In Sec. II, we introduce the effective chiral Lagrangians up to 0(p 3 ). In Sec. Ill, 
we calculate the EM decay amplitude and width order by order. In Sec. IV, we determine the low energy constants 
(LECs) using the quark model. In Sec. V, we discuss the numerical results and summarize. We collect some useful 
formulae and the numerical results of the EM decay of the bottom baryons in the appendix. 


II. THE CHIRAL LAGRANGIANS 


A. The Leading Order Chiral Lagrangians 


The three light quarks u,d,s form the fundamental representation under the flavor SU(3) transformation. The 
heavy baryons can be classified into multiplets in flavor SU(3) space. Now we concentrate on the ground states of the 
charmed baryon. When the two light quarks are in the flavor antitriplet, the quantum number of the heavy baryon is 
J p = i + . When the two light quarks are in the symmetric flavor sextet, the quantum number of the heavy baryon 
can be either J p = | + or | + . 

In this work, we investigate the heavy baryon systems employing the Chiral Perturbation Theory (ChPT) [22| . 
which provides a systematic expansion of the physical observables in terms of small momentum p and the mass of 
Goldstone bosons m. ChPT has been widely used to study the low energy hadron interaction. In the early stage, it 
was employed to study the purely mesonic system [23|, HH . Later it was extended to discuss the baryon-meson system 
[25l - [28l ] . In order to deal with the heavy baryon system, the heavy baryon chiral perturbation theory (HBChPT) 
was developed [32l-[35j , which provides a convenient framework to make a dual expansion in terms of both the small 
momentum and where M is the heavy baryon mass. 

In the flavor space, the pseudoscalar meson fields and spin-4 baryon multiplets are defined as follows 


f-° + js* 

■\Z2tt 

v V2K~ 


r' 


-7T° + -U 
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The spin-1 baryons B are the so-called Rarita-Schwinger vector-spinor fields [36], which are similar to Bq. In the 
framework of ChPT, one can construct the Lagrangians order by order while the specific building blocks satisfying 
the chiral transformation are defined as follows [29|, |3Ci[ . 


U(x) = e^ Hx \ u 2 = U 

D„B = d^B + T^B + BTl 
V m J7 = d^U-ir^U + iUl^ 

l\ \ 0 \ 

V = ^ = -eQAp, <9= 0-3 0 

\0 O' -\) 

The superscript T denotes the transpose in the flavor space. The pion decay constant F 0 
meson Lagrangian in the lowest order 0{p 2 ) is 

£( 2 ) = ^TV[V M [/(V M C/)t] 
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The leading order pseudoscalar meson and heavy baryon Lagrangians at 0{p 1 ) read [13l 

4 1} = ^Tr[B- 3 (Hp-M- 3 )B- 3 )+Tr[B 6 (0-M 6 )B 6 ] 

+Tr {B^l-g^i# - M*) + + 7 V D„) - 7 „(*# + M S)l u]B?} (5) 

Ant = SiT^A^s-Be) + 92 [Tr(i?6^75^3) + h- c -] + 53 [Tr(-B 6 M u M -B 6 ) + h.c.] 

+g4Tr(B^ Ufl B- 3 ) + h.c.] + g 5 Tr(B* 6 y l5 B* 6l/ ) + g 6 TiiB 3 ^ 5 B 3 ) ( 6 ) 

where C q 1 ^ and Ant are the free and interaction part. From the quark model and flavor SU(3) symmetry, the 

axial coupling constants ge = 0, gi = 0.98, 31 = — ^/§32 0]- The heavy quark spin flavor symmetry leads to 

the following relations among these coupling constants, i.e., g 3 = -^ 31,35 = —§ 31,34 = —V^> 92 - Numerically, 
31 = 0.98,3 2 = - 0 . 60,35 = -1.47 ,33 = 0 . 85,34 = 1.04. 

In the framework of HBChPT, the baryon field B is decomposed into the large component M and the small 
component 'H. Their relationships are 

B = e~ iMv - x {N + H) 

A (■ = e iMvx\+J_ B , n = e iMvx}_Z_l B ^ 

where v M is the on shell velocity. For the spin-§ baryon, the large component is denoted as T M . Such a decomposition 
is equivalent to separating the off-shell momentum k from the whole momentum p of the heavy baryon while the on 
shell part is Mv. The reduced effective Lagrangian only relies on the large component. Omitting some details, the 
nonrelativistic Lagrangian in the heavy quark limit reads 

£(!) = iTr[A r 3 {iv ■ D)U 3 ] + Tt[M 6 (iv ■ D - c^Ms] + Tr {f p [-g pa {iv ■ D - 5 2 )}T a } 

+2giTr(J\f e S ■ uJ\f 6 ) + 2g 2 [Ti(J\f 6 S • uj\f 3 ) + h.c.] + 3 3 [Tr(f l u u,_ l J\f 6 ) + h.c.] 

+ 34 [Tr) + h.c.] + 2g 6 Tt(T v S • u%) + 2g 6 Tr(JV 3 S ■ uM 5 ) (7) 

The mass difference parameters are defined as d-\ = Mq — M 3 ,S 2 = Mq* — Mq. In the isospin symmetry limit, 
<5i = 126.52 MeV, 5 2 = 67.03 MeV. We have not kept the jj corrections to the Lagrangian because their contributions 
are of higher order. 


B. The higher order Lagrangians for the EM decays 


Considering the conservation of the electric charge and energy, there exist twelve electromagnetic decay channels, 
which can be classified into three types: 

B 6 B 3 + 7 : 


E+ A+ + 




77/0 


C+ 


"0 


7 

-7 

7 


b; -> B 6 + 7 : 


e:++ e++ + 

e;+ 

-+ Z + + 7 

s:° 

^+7 

'—c 

' -»• S' c + + 7 

-^*/0 

'—c 

5' c ° + 7 


^B°+7 





4 


Bq —>■ B 3 + 7 : 


^*/o 


^ A+ + 7 
-7 
7 




2*' + ->. S + 


rjO 


The leading order Lagrangians in Eqs. ©-(O do not contribute to the above EM decays. So the lowest order EM 
decay diagrams arise from the 0(p 2 ) Lagrangians at the tree level. If we consider the next leading order corrections, 
we need the contribution from both the tree-level 0(p 3 ) Lagrangians and the chiral loop corrections from the leading 
order Lagrangians. 


The higher ordered Lagrangian can be constructed with the help of the new building blocks f^ v . 
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for the charmed baryons, Qb = 



for the bottomed baryons 

( 11 ) 

The transformation property of the building blocks can be seen in Ref. [3l| . What we need in the following is that 
the operator f±, transforms as the adjoint representation. We notice that the direct product of 3 x 6 = 8 +10 contains 


only one adjoint representation. Therefore, there is one independent interaction term in the 0(p 2 ) Lagrangians for 
the EM transitions between the 3 and 6 / 6 *-dimensional baryons. 

In contrast, the direct product of 6 x 6 contains one singlet and one octet. Accordingly, there are two independent 
interaction terms for the EM transitions between the 6 and 6 *-dimensional baryons. The octet piece corresponds to 
the /+, term while the singlet piece corresponds to the trace of /+ in Eq. (fl2l) . The 0(p 2 ) Lagrangians are as follows 


'4J = ATrR^V+^ + h.c. 

< 4 2 6* = ^-Tr B 6 e^f+iD p B 6 ; + A- n B 6 ^(/+ )iD p B 6 J + h.c. (12) 

46 * = ^TrB- 3 e^f+„iD p B 6 ; + h.c. 

where Tr and () denote the trace in the flavor space, and A x = AttFo. The coefficients fi s are the new dimensionless 
LECs. The 8 -dimensional operator /+, and the 3, 6 -dimensional baryon fields couple into the possible flavor singlet. 
Similarly, the 0(p 3 ) Lagrangians for the EM decays are constructed as 


/•(&) _ 

^66* — 


||TVR 6 V A /+7 A 75^ [M i?6i + ^TrBeVxf+^iD^B*. + ||tVR 6 V a /+ iD x ^ l5 B”l + h.c. 


(13) 


+ ^TrR 6 V A (/+) 7 A 7 5 zD [^ 6 1 + ^TrB 6 V x (f+)^ 5 iD^B%, + ||t^R 6 V a (/+ + h.c. 


46* = ^TrR3V A /4 7 A 75^[^ 6 1 + ^TrR3V A /+ !/ 7 [M 75^ 1 B 6 A . + ||TrR3V A /4* J D A 7^7 5 R 6 1+h.c. (14) 


where D^B U ^ = D^B 1 ' — D^B^ and the other commutators are similar. The coefficients hi s are also dimensionless 
LECs. 

Because of the constraint of the above flavor structure, we notice that there does not exist EM transitions between 
the neutral 3 and 6 / 6 *-dimensional baryons at the tree level up to 0{p 3 ). This seemingly surprising result is the 
manifestation of the famous flavor SU(3) U-spin symmetry, which interchanges the strange and down quark. Notice 
that the s and d quarks carry the same charge. With the exact SU(3) flavor symmetry, the s and d quarks are identical 
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particles both in QCD and QED! Therefore, the Pauli principle forbids the s and d pair within the 2° baryon to form 
the flavor antitriplet. In other words, the decay processes 2(? — > 2° + 7 and 2*' —► 2° + 7 are forbidden at the tree 
level up to 0(p 3 ). 

If we consider the explicit SU(3) symmetry breaking, we can construct the tree-level Lagrangians at 0(p 4 ) such as 

£ (4) cx TrB~ 3X+ a^(f+)B 6 + h.c. (15) 

where x+ is the building block relevant to the current quark mass [3l|. Since we truncate the chiral expansion at 
0(p 3 ), we do not try to exhaust all the possible 0(p 4 ) Lagrangians. 

On the other hand, the decay processes 2(° —► 2° + 7 and 2*' —> 2° + 7 can occur through the chiral loops at 
0(p 3 ) if we consider the explicit SU(3) symmetry breaking. Its divergence would be absorbed by the above tree-level 
Lagrangians at 0(p 4 ). We will include these corrections carefully in the following section. 


III. THE EM DECAY WIDTHS IN HBCHPT 
A. The EM Decay Widths at 0(p 2 ) 


We calculate the EM decay amplitude in Coulomb gauge. For a nonpolarized process, the summation of the photon 
polarizationyields the metric h^ v 

E r* r _ 7 

€ li e v ~ —/ V*' 

r= 1,2 

We adopt the metric h^ v in Ref. [38j |. 

k l * “1“ Iflll/ 

= 9 ^ (/. n) 2 - P 

where n = (1,0, 0,0) and l is the photon 4-momentum. We plot the lowest order Feynman diagrams from Ea. (|12l) in 
Fig. [U The single line represents a spin-i baryon and double line a spin-| baryon. The solid dot represents an 0(p 2 ) 



FIG. 1: The 0(p 2 ) tree level process 


EM vertex. 

Summing the photon polarization r, the spin directions of the initial and final baryon states Sa, Sb , we obtain the 
square of EM decay amplitude modular for the spin-i to | and spin-| to | transitions respectively. 
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(17) 


where M A ,M B represent the baryon mass of the initial and the final states, t is the modular of the photon 3- 
momentum. In the rest frame of the initial state, t only depends on the baryon mass. 

( M a + M b )(M b -M a ) 

2 M b [ ’ 

For each channel, the coefficients ’s are listed in the 2nd column of Table Q] 

Since some LECs always appear in a combined form in the decay amplitudes, we abbreviate hij = hi + hj where 
hi etc arises from the tree-level Lagrangians at 0{p 3 ). We also collect the LECs h r . } = /q + hj in Table Q] For the 
bottomed baryon systems, the LECs / b , h b are marked with the superscript b for distinction. The corresponding 
coefficients are collected in Table HT1 
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TABLE I: The coefficients D^\ D^\, D^\ in the EM decay amplitudes. 

EM decay processes Coeff.P^ of Q(p 2 ) CoeffiD^ of Q(p 3 ) CoeffiD^ of Q(p 3 ) 
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TABLE II: The coefficients D ( ^, , 

/o') 

D\j 2 for the bottomed baryons. 

EM decay processes Coeff.D+ of 0(p 2 ) Coeff.D ) -+ ) 1 of 0(p 3 ) Coeff -D\ 3 ^ 2 of 0(p 3 ) 
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B. The 0(p 3 ) contributions at the tree level 

The next leading order contributions arise from the tree-level Lagrangians in Eas. (1131141) and chiral loop corrections 
from Eqs. ©-©• The Lagrangians in Eos. (1121131) contribute to the tree-level process between the spin-^ and | baryons 
in Fig. [2] In fact, there are two other Lagrangians involving /+, at 0(p 3 ). One is Bj^V 1 ' f pu B which leads to the 
EM decay amplitude with the factor — Z 2 e (U . Under the Coulomb gauge, this term vanishes. The other possible 

interaction Lagrangian is B^f^D p B a . In the leading order of the heavy baryon expansion, its contribution also 
vanishes since S!^T„ = 0. 

There are two Lorentz structures from Fig. [2] To sum up, all the EM decay amplitudes at the tree level up to 
0{p 3 ) can be written as 

| iM ( 2) = D^Ms(l^ v - U^e^VpU^uf* 

\ iM (3) = w[ 3) Mbv ■ lu SA S ll 'l lt iea’uf? + iD {3) M B v ■ Iu Sa S^' e„4 a >u^ 
where the Pauli-Lubanski spin operator S M is defined as S^ = — 575 ( 7 + — v M ). The coefficients D^\ D^’s are also 
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FIG. 2: The 0(p 3 ) tree level process. 


listed in Table Q] In Coulomb gauge, the cross term between the 0(p 2 ) and 0(p 3 ) EM decay amplitudes reads 

E ( iM ( 2 ) )*iM (3) + (iM { 3 ) )*iM {2) = -y MaM^ 3 D&(D^ - D < 3) ) (20) 

r,S A ,S B 


The square of the decay amplitude from the 0(p 3 ) diagram is 


E (l?Xi| 2 + \iM 2 \ 2 ) 

r,S A ,S B 




M A M 3 B t 


( 21 ) 


C. The 0(p 3 ) chiral loop contributions 



FIG. 3: 0(p 3 ) loop level process 

At the one-loop level, there are six possible Feynman diagrams as shown in Fig. [3] All the vertices are from Eqs. 
(1415161) . In diagram a, the photon vertex is from the tree-level EM interaction term while the meson vertex from the 
strong interaction terms. In diagram b, the photon-meson vertex is from the chiral connection. In diagram c and d, 
both two vertices are from the strong interaction and seagull terms. In diagram e, the meson vertex arises from the 
chiral connection. 

The external baryons in diagrams b and e must be the same. So they do not contribute to the EM decay process. 
Diagrams c and d vanish in the heavy baryon mass limit. The nontrivial recoil corrections appear at 0(p 4 ) which is 
beyond our chiral truncation order in this work. The decay amplitude from diagram a between spin ^ to ^ baryons 
is proportional to either S ■ v = 0 or v ■ T = 0. Hence its contribution vanishes. Similarly the decay amplitude 
between spin i to \ baryons from diagram a is always proportional to e • v. Under Coulomb gauge, this amplitude 
also vanishes. 

Diagram f can be classified into two different types according to the Lorentz structure in the loop integrals, which 
are displayed in Fig. 0J For g$ = 0, the contributions from two of these Feynman diagrams in Fig. |1] vanish. 



















FIG. 4: The flavor structures of loop-f. 

1 . spinto ^ diagram 

For the 6 —> 3 EM decays, the decay amplitude from the chiral loops can be written as 

iM lo ° p ' f = iM\ + iMu = e^ a0 v a u s i A Spu 83 - l^e^) l -b\- \ibn 

J 2 6 


The function [|&i — |i 6 u] contains all the coefficients involving gi etc and the loop integrals. The loop-f Lorentz 
structure is the same as that of the 0(p 2 ) tree diagram. The divergence of the chiral loops can be absorbed by the 
0(p 2 ) LECs. We list the functions of the three channels below. 
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The baryon mass splitting and the meson type in the loop can be read from the functions directly. The function 
7 i(m, 6 , t) arises from the loop integral, which are defined in Appendix IVIII 
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2 . spin -1 to i diagram 

For the spin-| to | EM decays, the Lorentz structure of the loop-f amplitude is different from the structure of the 
0(p 2 ) tree level one, but is the same as that of the 0(p 3 ) tree level one. Both type-I and II loop-f amplitudes have 
two structures, which are marked with the subscript 1 , 2 . 
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The function [2b\ + 2bu) contains all the coefficients involving the coupling constants and the loop integrals, too. 
For the six EM decays between the 6 and 6 *-dimensional baryons, we have 
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F 2 

r o 


F 2 

r o 


H 3 ( m K , 0 ,£) | *7s (m^,0 ,£) + *73 (m^,0 ,£) | *7 8 (m w ,0,£) 
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oi . , I _ l 'i'Y3{rn K ,S 2 ,e) , i'ya (m K ,S 2 ,£)\ , 

26n + 26m| n o < _ n ,o - 9i93e ^^- 1 - ^2 -j+9294e^ 


/ 2z73 (m Kl 5 3 ,£) 2 z 7 9 (m K ,S 3 ,£) 


+9395 e 


0 x 0 
/ *73 (tox, 0 , £) 2Z73 [ m K , 0, £) 




3+o 


3+o 


2&I2 + 2frn2 | n o^o.o = 
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F 2 
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For the three EM decays between the 3 and 6 *-dimensional baryons, we have 

il 3 (rnK,S 2 ,£) i^ g (m K ,S 2l £) 2iV2j 3 (m w , 5 2 , l) 2iV2 r y g {m n ,S 2 ,£) 
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+n 
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3 v^+o 2 


+ 
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3+o 
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v^+o 2 
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+9495 e 

26ii + 26in|-+ -*,+ = 92936 — 

0 0 V 

+9495 e 


*73 (*Ffc:,0,£) *7s(TOif,0,£) 2z v / 273 (m T , 0, £) 2iy/2'y s (m*-, 0, £) 


v^+ 0 2 

3V2+0 2 

F 2 

r o 

3+o ) 

2*v / 273 ( niK,S 2 ,£) 

2iy/2^ g (rriK,6 2 ,£) 

*73 (m n ,S 2 ,£) 

*79 (m n ,S 2 ,£) 

F 2 
r o 

F 2 
r o 

V2F 2 

V2F* 

4iy/2j 3 (m K , 0,£) 

2iV2'y 8 (mK,0,£) , 

iV 273 {m„, 0,1) 

*78 (m n ,0,£) 

3+o 

3+o 

00 

3V2+o 2 


2&12 + 26112| H +, H ./+ = 9293e ( — 


2iV2'y 8 (m K ,5 2 ,£) ij 8 (m v ,d 2 ,£) 


F 2 

r o 


v^+ 0 2 


+9495 e - 


2 z \/273 (tok,0,£) 2iV2'y 8 (m K ,0,£) * 7 3 (m,r, 0 ,£) *7s (m w , 0,£) 


+o 2 


3+o 


V2+ 0 2 


Sv^+q 


26 n + 26 ih| 30< _ s ./o = 92936 ( — 


*73 (m K , S 2 ,£) _ *79 {m K ,6 2 ,£) *73 (m„,S 2 ,£) *79 {rn^,5 2 ,£) 
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iV2j 3 (m K ,0,£) * 7 s(mx, 0 ,£) *\/273 (m w , 0, £) *78 (m w ,0,£) 


26 i 2 + 26 ii 2 | HO< _ s ./o = 92936 


+9495 e 

/ ij 8 (m„,d 2 ,£) i'y s (m K ,6 2 ,£) 


3+o 


3-\/2+o 


3+o 


3-\/2+q 


V v^+ 0 2 v^+o 2 


, *73 (*Ar:,0,£) *7 8 (m-if, 0,£) *73 (m„., 0, £) *78 (***•*■, 0, £) 

+9495e-^ 33 ;- . -1- ^333 -h 


v^+n 2 


3-\/2+o 


v^+o 2 


Sv^+o 2 


IV. THE DETERMINATION OF THE LOW ENERGY CONSTNATS (LECS) 


In principle, the low energy constants (LECs) in the chiral Lagrangians should be extracted through fitting to the 
experimental data. However, the available information on the heavy baryon system is still rather poor, which renders 
the extraction of the low energy constants /, h nearly impossible. We try to estimate the LECs using the quark model. 
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At the quark level, the flavor wave functions of the spin i baryons can be written as follows. 


A °~ = du ) c 
Sj" = —— (us — su)c 

V2 

s c = ^=( ds-sd)c 


E+ + = uuc 

E+ = —-=(ud + du)c 
v 2 

E° = ddc 

E!+ = 4 ={us + su)c 

S'° = -^={ds + sd)c 
. = ssc 


( 22 ) 


The flavor wave functions of the spin-^ baryons are similar. 
The spin wave functions of the heavy baryons are 

1 


= yf(xix-ix-| - x-ixix-i) 

X|| = vfxixix-i - ^(XIX-IXI +X-iXiXi) 

1 , . 

= y|(XiX-|X-| + X-iXiX-i) - y gX-iX-iXi 

Xfl =X 1 X|X| 
x#-| = X-iX-iX-i 


(23) 


Comparing the matrix element of the same EM decay process at both the hadron and quark level, one can express 
the LECs in terms of the constituent quark masses and electric charge. For example, let’s consider the E+ —> 7 + A+ 
decay process. At the hadron level, the Lagrangian is 


/•( 2 ) 

' C E+-> 7 +A+ “ 


2V2 

AT 


/ 2 eA+a^E+F M1/ 


(24) 


In the rest frame of the initial state, we choose the momentum of photon along the 3rd axis with the polarization 
along the y-axis. The matrix element of such a decay reads 


A+ —- 


A, 


4V2 

^7 


S 2 & 


(25) 


where we have dropped the photon momentum and the image imaginary unit etc. These common factors also appear 
in the matrix element at the quark level. On the other hand, we consider the same matrix element at the quark level. 
The Gordon decomposition separates the QED interaction Lagrangian into the electric and magnetic parts. The light 
quark magnetic interaction reads 


Am = --A 77 —— wcAA - - -FscAa) (26) 

4 \ om u Srrid orn s J 

For the heavy quarks we have 

Am, c = —e-^—ca^cF^, Am,b = ^e-^—ba^bF^ 

4 oTflc 4 0777/5 


respectively. We notice that the heavy quark magnetic interaction only contributes to the EM transitions between 
the 6 * and 6 -dimensional baryons. 
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The flavor-spin wave functions of A+ and X+ at the quark level are 


A+,-^y = -j={\udc) - \duc)) ® 
= ^=(Mc) + \ duc » 5 


1 

L 7 I 

2 


(I m> -1 m» 


3' m >-V6 (,nt) 


Itt)) 


(27) 


The magnetic interaction and tensor operator a flips the spin of the up and down diquark system. The matrix 
element at the quark level reads 



x+ 1 

Lc ’ 2 


— ( 1 i M 

3x/3 \m u 2 m d ) 


With the commonly used values of the constituent quark masses 


m u = 335.5MeV, mj = 339.5MeV, m s = 510MeV, m c = 1550MeV, mb = 4730MeV, 


one obtains the value of the LEC fa 


(28) 


fa = 0.166 

The extraction of the other LECs for the 6 * —> 6 process is more complicated. Let’s take the decay channel 
S*+ X+ + 7 for an example. Both the 0(p 2 ) and 0(p 3 ) Lagrangians will contribute to the matrix element at the 
hadron level. 


£ ( 2 ) 


-(* + 


f 3 - 7kT^) /3 eM s , + 


V2A 2 X V2A 2 X 


S+ I 
Lc, 2 


£ 0 ) 


sr. |) = ) e(« a+ + M 2 p ( M 2 .. - M 2f ) 


The matrix element at the quark level reads 


x+ 1 

C ) 2 


c 

y*+ 3 \ 

- v 

f 1 

1 1 ) 


c ’ 2 / 

3x/6 

\m u 

2rrid m c ) 


(29) 


There are similar equations at both the hadron and quark level for all the other 6 * —> 6 processes. 

At 0(p 2 ), there are only two unknown LECs fa and fa while there are more equations. The best fit in the SU(3) 
flavor symmetry limit leads to 


fa = 0.099, fa = -0.037. 

In the extraction of the 0(p 3 ) LECs h \3 etc, we use the above values as input and take into account the explicit SU(3) 
symmetry breaking, i.e., using the experimental values of the heavy baryon masses and different values for m u and 
m s . Solving these above equations simultaneously, we get the values of the LECs at 0{p 3 ) 

h 13 = 0.073, h 23 = 0.034, h 13 = 0.035, h 23 = 0.061 

Similarly, we can extract the LECs for the 6 * —> 3 process. 

fa = 0.133, /i 46 = 0.056, h 5 e = -0.039 

All the values of LECs are collected in Table m 
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TABLE III: The values of the dimensionless LECs for the charmed baryon systems from the quark model. 



0(p z ) tree level 0(p* 

) tree level 

h 

0.166 

X 

h 

0.099 

X 

h 

-0.037 

X 

U 

0.133 

X 

hi3 

X 

0.073 

fl23 

X 

0.034 

hl3 

X 

0.035 

h 23 

X 

0.061 

^46 

X 

0.056 

^56 

X 

-0.039 


TABLE IV: The EM decay width of the charmed baryons in unit of keV. 


EM decay process 0(p 4 ) tree level 

0(p 5 ) cross 

term 0(p 6 ) tree level 

0(p b ) Cross term with 0(p°) Tree+loop 

loop contribution correction 

Full width 

S+ 7 + A+ 

149.41 

X 

X 

X 

164.16 

164.16 

3c + —► 7 + H 4 

43.26 

X 

X 

X 

54.31 

54.31 

S' c ° —> 7 + S° 

X 

X 

X 

X 

0.02 

0.02 

S* ++ -5-7 + SJ+ 

8.21 

1.46 

0.13 

2.88 

0.51 

11.60 

S*+ -»• 7 + E+ 

0.35 

0.23 

0.08 

0.32 

0.18 

0.85 

E*° —>■ 7 + E° c 

2.93 

-0.41 

0.03 

-0.03 

0.02 

2.92 

sr + -> 7 + H'+ 

0.48 

0.34 

0.13 

0.43 

0.19 

1.10 

K'° -»• 7+s; 0 

3.87 

-0.58 

0.05 

-0.04 

- 

3.83 

fi*° — > 7 + 

4.82 

-0.75 

0.06 

-0.03 

0.03 

4.82 

S* + -»• 7 + A+ 

756.40 

18.43 

7.23 

125.74 

10.83 

892.97 

H* ,+ ->7 + 2+ 

402.92 

7.63 

2.33 

94.70 

4.49 

502.11 

2^'° ^7 + 2° 

X 

X 

X 

X 

0.36 

0.36 


V. NUMERICAL RESULTS AND DISCUSSIONS 


So far, all LECs and the chiral loop corrections have been obtained. The EM decay width reads 


r = 


i 

8irM\ 


\M \ 2 


(30) 


We collect the numerical results of the EM decay widths of the ground state charmed baryons in Table HVl (in unit of 
keV). The notation ” x” means no contribution in this channel and means the contribution is strongly suppressed. 
In Table ITVl the first three rows, middle six rows and the last three rows correspond to the 6 —> 3, 6 * — > 6 and 6 * — > 3 
EM decays respectively. 

The EM decay width is proportional to the cubic of the photon momentum, which depends on the mass difference 
of the initial and final heavy baryons. The larger the mass difference, the larger the decay width. Therefore the 
6 */ 6 —>• 3 EM decay widths may reach several hundred keV while the 6 * —> 6 transition width is only a few keV. 

From Table IIV1 the numerical results of the EM decay widths show very good convergence of the chiral expansion. 
For the 6 — > 3 decays, the leading order widths are much larger than the loop contributions. For the nine spin-| to 
t transitions, the leading order amplitude yields the EM decay width at 0(p 4 ). The cross term from the tree-level 
0(p 2 ) and 0(p 3 ) amplitudes contributes to the decay width at 0(p 5 ) while the 0(p 3 ) amplitude alone contributes 
to the widths at 0(p 6 ). From the second, third and fourth column of Table IWl the numerical results of the 0(p 4 ), 
0(p 5 ) and 0(p 6 ) corrections decrease accordingly, consistent with the chiral power. When the loop contributions are 
added in the fifth and sixth columns, the numerical values also maintain the same decreasing trend. 

A remarkable feature in Table m is the strong suppression of the 3(? —>• 7 + and E *' —> 7 + S° decays. As 
we have pointed out before, the exact SU(3) U-spin flavor symmetry forbids the two neutral EM decay processes 
S(.° -t 7 + “S and S * /0 —> 7 + 3°. There do not exist tree-level contributions at 0(p 2 ) and 0(p 3 ). The non-zero 
tree-level contribution to the decay amplitude starts at 0(p 4 ) which explicitly breaks the SU(3) flavor symmetry 
through the current quark mass matrix. On the other hand, the chiral loop corrections appear at 0(p 3 ) if we use 
the physical mass for the eight Goldstone bosons. In fact, the chiral corrections from the pion, kaon and eta mesons 
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will cancel exactly at 0(p 3 ) if their masses and decay constants are the same, which is also required by the U-spin 
symmetry. From Table IIVI we do notice that the numerical values of the loop contributions of these two channels are 
indeed much smaller than their counterparts in the same multiplet. 

In short summary, we have investigated the electromagnetic decays of the heavy baryon systems within the heavy 
baryon chiral perturbation theory. After introducing the leading order chiral Lagrangians and constructing the 0(p 2 ) 
and 0(p 3 ) Lagrangians, we have systematically calculated the twelve EM decay amplitudes up to 0(p 3 ). With the 
help of the quark model, we estimate all the low energy constants at 0(p 2 ) and 0(p 3 ). The convergence of the chiral 
expansion is quite good. 

In the future, the EM decay widths of the charmed and bottomed baryons may be measured at facilities such as 
LHCb and JPARC. We sincerely hope our present work will be useful to future experiments. Moreover, the explicit 
chiral structures of the heavy baryon decay amplitudes derived in this work may also be useful to the possible chiral 
extrapolations of the future lattice simulations of these EM decay amplitudes. 
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Appendix 

VI. EM DECAY WIDTHS OF THE BOTTOMED BARYONS 

In the same framework, we estimate the EM decay widths of the bottomed baryons. The chiral Lagrangians of the 
bottomed baryons are the same as those for the charmed baryons. All the coefficients f b , h b are collected in Table HT1 
They can also be determined from the quark model. The value of the LECs are listed in Table [Vj 


TABLE V: The dimensionless LECs for the bottomed baryon systems. 


0(p z ) tree level 0{p 5 ) tree level 

n 

0.166 

X 

fa 

0.047 

X 

fa 

-0.021 

X 

ft 

0.058 

X 

hi3 

X 

0.369 

h-23 

X 

-0.058 

hi3 

X 

-0.096 

h-23 

X 

0.068 

/l46 

X 

0.087 

hie 

X 

-0.061 


At present, the experimental data on the bottomed baryons is scarce. We choose the average value of the bottomed 
baryon mass as M| = 5891.87MeV Q. We naively assume that the the mass splitting between 6 and 6 *-dimensional 
is inversely proportional to the heavy quark mass. 

(■ M b , - M b )m b = (M 6 c , - M 6 c )m c 
Then we obtain the mass splitting for the bottomed baryons as 

= 21.97MeV 

We let 8 b ,8 b take roughly the same values as in the charmed baryon systems since these splittings arise from the 
strong interaction of the light diquark system. The above values will be used in the calculation of the chiral loop 
corrections. 

Since the EM decay widths are very sensitive to the photon momentum, we explicitly keep the photon momentum 
t in the decay width formula which shall be measured by future experiments. The EM decay widths of the bottomed 
baryons are listed below in unit of keV while the photon momentum £ is in unit of MeV. 


r sg 

0-0 

< 

+ 

(- 

T 

= 4.15 

X 

nr 







S* 

[I] 

f—1 

^ 7 +Hg = 

= 4.67 

X 

10" 

-H 3 






r H /- 

^b 

-^7+3“ 

= 2.18 

X 

10" 

-8 £ 3 






r £* + ^ 7 +£+ = 

= 5.37 

X 

10" 


+ 

5.74 x 10" 

7(4 

+ 3.20 x 

10“ 9 £ 5 

r y]*o 

^b 

^7+Sg = 

= 5.12 

X 

10" 

-6 £ 3 

+ 

3.21 x 10" 

8(4 

+ 2.89 x 

10 -W(5 

r E *-_ 

^b 

-^7+S^ 

= 7.86 

X 

10" 

-6(3 

+ 

13.10 x 10 

-8( 

4 + 14.22 

x 10" 10 £ 5 


= 5.30 

X 

10" 

-6(3 

+ 

2.98 x 10" 

8(4 

+ 2.14 x 

1 0 " 10 ^5 

r„ , = 

i E*'-^ 7 + S' b 

= 7.86 

X 

10" 

-6(3 

+ 

13.82 x 10 

-8( 

4 + 13.64 

x 10~ 10 £ 5 

r n*- 

^■7+n 6 

= 8.29 

X 

10" 

-6(3 

+ 

15.59 x 10 

-8( 

4 + 26.86 

x 10" 10 £ 5 

r s;° 

->7+Ag = 

= 7.55 

X 

10" 

-5(3 

- 

1.85 x 10" 

7(4 

+ 2.31 x 

10“ lo £ 5 

r^,o 

l—, b 

^7+3« = 

= 7.95 

X 

10" 

-5(3 

- 

3.13 x 10" 

7(4 

+ 2.84 x 

1 0 “ 1o £5 

r_»,_ 

“b 

-^7+3“ 

= 2.15 

X 

10" 

-11(1 







We also list the numerical results of the EM decay widths in Table fVTl where the masses of the initial and final state 
bottom baryons are known experimentally pj. 
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TABLE VI: The EM decay widths of the bottomed baryons decay in unit of keV. 


EM decay process 0(p 4 ) tree level 

O(p 0 ) Cross term with 0(p b ) Tree+loop 

loop contribution correction 

Full width 

ES -> 7 + Ag 

255.90 

X 

287.65 

287.65 

^7+s 6 + 

0.48 

0.11 

0.01 

0.60 

Ef^ 7 + E“ 

0.04 

0.01 

- 

0.05 

7 + 

0.06 

0.02 

- 

0.08 

sr “ 7 + 

701.63 

-361.29 

94.83 

435.17 


281.32 

-168.70 

23.34 

135.96 

E b'~ -!-7 + S^ 

X 

X 

1.87 

1.87 


VII. INTEGRALS AND FUNCTIONS 


1. The integral with one meson line 


7(m) = 


167T 2 


R + In ■ 


2. The integral with one meson line and one baryon line 


J(m, uj) 


1 

8 ^ 


uj(R — 1) + ui In — tt + K 
l-i 


, ui = v ■ k + S 


(31) 


and 


2\J u: 2 — m 2 arccosh- 2 i 7 n/ uj 2 — to 2 , uj > m 

m 

K = < — 2 \j uj 2 — m 2 arccosh-, uj < —to 


2\/ to 2 - uj 2 arccos-, uj 2 < m 2 

m 


3. The integral with two meson lines and one baryon line 
3.1 scaler integral 


L = 


1 f dr q 


i J (2n) d ( m 2 — (q + l) 2 ){m 2 — q 2 )[v ■ q + uj\ 


l 2 = 0 


(32) 


L = 


1 


8tt 2 v ■ l 


farccosh-—-^ farccosh- 


- 7 xi 


In ( \/(lo — v ■ l) 2 — m 2 + uj — v ■ l ) — In ^\/ uj 2 — m 2 + uj ] 

, /3 < — to 


')]’ * 


> TO 


( arccosh-—- 1 — (arccosh- 


( —uj + v ■ l\ ( —uj 

arccos- — arccos- 

\ to ) \ TO 


(33) 


/3 2 < to 2 


where/3 = uj — v ■ l 
3.2 vector integral 




1 

i 


d d q _ q ^ _ 

( 2ir) d (m 2 — (q + l) 2 )(m 2 — q 2 )[v ■ q + uj\ 


7 i^ + 72 


(34) 


= - -{him) — wL{rn,w,v • l)) 

v ■ L 


1 


[J(to,/3) - J[m,w)\ 


7 i {m,u;,v ■ l) 


2 (vl) 2 


(35) 
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72 (m,w,v-l) = 


2 (vl)' 


;[J{m,0) - J(m,ui)} 


in the rest frame, v • l = t, which is the modular of photon 3-momentum. 
3.3 Second-order integral 


1 f d d q _ q ^_ 

i J (2n) d (m 2 — (q + l) 2 )(m 2 — q 2 )[v ■ q + u] 
= 73 9)1V + 74 lpl v + 75 Vplv + 76 l^V v + ^V ll V v 


7 3 (to,w,7) = - 

74 (m, w, 7) = 


1 


2(d — 2)£ 


{wJ(m, ui) + (j3 — 2u> + 2 £)J(m, j3) — 2m 2 £L(m, to, £)} 


1 


r{(w — dui)J(m , w) + ./(to, f3)(—(d — 1)(/? — 2w) — 2/) 


2(d-2)£ 3 

+2/ ((d — 2)w 2 + m 2 ) L(m , w, £) — (d — 2)£/ 2 (m)(2o; + /)} 


75,6 (m,w,l) = 


1 


2(d-2)£- 


-{ J(to, (3)(/3(d — 1) — dco + 21) + uJ(m , w) — 2 m 2 £L(m, u/, /)} 


1 


7 7 (m,u,£) = -{wJ(m,w)-/3J(ra,i3)} 

d arises from the dimensional regularization in d-dimensional space. 
3.4 Third-order integral 


(36) 


(37) 

(38) 

(39) 

(40) 

(41) 


r = I /'_A_ g/l gog /3 _ 

pa/3 ij (2ir) d (m 2 — (q + l) 2 )(m 2 — q 2 )[v ■ q + ui] 

= Isd^odp + 79g«3^a + hodden 

+lll9 l ia.Vp + 7l2g/y3U a + 7l3ga/3?V 

~t“7l4^/dod/3 4“ 7l5^/^o^/3 T 'YlQlfj.Valp + 7l7^/i^a^/3 

+7l8^aW/3 + llQV^laVp + "f20V^,V a lp + 721U Q Vp 

For completeness, we list all the functions 7 * as follows, although only j 8 and 79 are used in the calculations of 
loop-f. Some functions may be useful in the calculation of the magnetic moment of heavy baryons. 


78 , 9 ,io(m,u/,/) 


711,12,13 (m,ui,£) 


7i 4 (m,o/,/) 


715,16,17 (m, a/, £) 


718,19,20 ( m,UJ,i) 


721 [m, u), £) 


2{d-2)\d-l)P { - ^ 2 + (d “ 2)m2) J(m ’^ (43) 

+J(m,u/) ((d — 2)m 2 + uj(—d£ + uo + £))+£ (2(d — l)m 2 (him) — Leu) + (d — 2 )I(m))} 

2(d-2)(d-l)£ {J{m ’ ((d " 2)m2 + ' 8(2(d - ^ - P( d ~ 2 ))) ( 44 ) 

+J(m,u/) (—(d — 2)m 2 + (d — 2)o/ 2 — 2(d — l)u//) — (d — 2)£I(m)} 

6(d-2)(d-l)d^ { ~ 3a!J(m ’ ^ ^ 2((d ” 2)d + 3 ) + ( d - 2 ) (3(d ^ l)^ 2 + 2m 2 )) (45) 

+3 dj(m, oj) (2(d — 2)m 2 + ui (((d — 2)d + 3 )cj — 3(d — l) 2 /)) 

—6(d — 1 )dLuj£ ((d — 2)u/ 2 + 3m 2 ) + (d — l)£h(m) (l2(2d — l)m 2 + (d — 2)d (6a/ 2 + 2/ 2 + 3uj£)) 

+6(d - 2)((d - l)d + l)//(m)} 

2(d-2)(d- 1)7 3 ^ 2 ^ “ 1 ^ W ' 2 ^ “ /2 ( m )) - ( d - 2 )7(m)) - ((d - 2)m 2 + a/ 2 ) J(m, a/)) 

+J(m, /3) (2/3 2 + (d - 2) ((d - 1)1? 2 + 2m 2 ) + 2 /3(d - 1 )£)} (46) 

2(d-2)\d-l)P {J(rn ’ + d - 2 ) “I + (d - 2)m 2 - (d - 2)a/ 2 ) (47) 

+J(m, p) (d(/3(d - 2) - 4(d - 1)£) - (d - 2)m 2 ) + (d - 2) 2 (-^)7(m)} 

2(d- l)A ^ 2d ~ m2 ^ J< - m ’ ,8 ' ) + ^ m2 _ da;2 ) J< - m ’ + 


(48) 
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